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Abstract. Traditional black- hole binary puncture initial data is conformally flat. 
This unphysical assumption is coupled with a lack of radiation signature from the 
binary's past life. As a result, waveforms extracted from evolutions of this data display 
an abrupt jump. In Kelly et al pQ, a new binary black- hole initial data with radiation 
contents derived in the post-Newtonian (PN) calculation was adapted to puncture 
evolutions in numerical relativity. This data satisfies the constraint equations to the 
2.5PN order, and contains a transverse-traceless "wavy" metric contribution, violating 
the standard assumption of conformal flatness. Although the evolution contained less 
spurious radiation, there were undesired features; the unphysical horizon mass loss and 
the large initial orbital eccentricity. Introducing a hybrid approach to the initial data 
evaluation, we significantly reduce these undesired features. 
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1. Introduction 

The field of Numerical Relativity (NR) has progressed at a remarkable pace since the 
breakthroughs of 2005 [21 El S] with the first successful fully non-linear dynamical 
numerical simulation of the inspiral, merger, and ringdown of an orbiting black-hole 
binary (BHB) system. In particular, the "moving-punctures" approach, developed 
independently by the NR groups at NASA/GSFC and at RIT, has now become the 
most widely used method in the field and has been successfully applied to evolve 
generic BHBs. This approach regularizes a singular term in the space-time metric 
and allows the black holes (BHs) to move across the computational domain. Since this 
breakthrough, BHB physics has rapidly matured into a critical tool for gravitational 
wave (GW) data analysis and astrophysics. Recent developments include: studies of 
the orbital dynamics of spinning BHBs El [7J El EE QUI ELT], calculations of recoil 
velocities from the merger of unequal-mass BHBs [121 H31 E] , the surprising discovery 
that very large recoils can be acquired by the remnant of the merger of two spinning 
BHs [151 HS1 HZl HH H91 EQl EU [221 [231 EH [251 E6J ETl EJ EH EHl EU], empirical models 
relating the final mass and spin of the remnant with the spins of the individual BHs 
[3lT 1321 1331 IM1 1351 1361 1371 138] , and comparisons of waveforms and orbital dynamics of 
BHB inspirals with post-Newtonian (PN) predictions [391 SHI SH S21 SHI S31 SSI SSI ST] . 

One of the important applications of NR is the generation of waveforms to assist 
GW astronomers in their search and analysis of GWs from the data collected by ground- 
based interferometers, such as LIGO [33] and VIRGO [33], and future missions, such 
as LCGT [50], LISA [51], ET [52] and DECIGO 153 . BHBs are particularly promising 
sources, with the final merger event producing a strong burst of GWs at a luminosity 
of Low ~ 10 22 lJ3, greater than the combined luminosity of all stars in the observable 
universe. The central goal of the field has been to develop the theoretical techniques, 
and perform the numerical simulations, needed to explore the highly dynamical regions 
and thus generate GW signals from a representative sample of the full BHB parameter 
space. Accurate waveforms are important to extract physical information about the 
binary system, such as the masses of the components, BH spins, and orientation. 

With this in mind, we note a drawback shared by most present-day comparable- 
mass black-hole binary simulations in three spatial dimensions: they are performed 
using conformally flat initial data given by the Bowen-York (BY) [51] prescription as 
applied by Brandt and Briigmann [55J. This prescription, while numerically convenient, 
lacks physical realism. The unique stationary vacuum black-hole solution to Einstein's 
equations is the Kerr solution, which is not conformally flat for non-zero spin. We 
cannot approximate this with single-puncture spinning BY data without also including 
unphysical radiation. 

Conversely, approximating an inspiralling black-hole binary system with two- 
puncture BY data will leave out the gravitational radiation expected in physical 

| This luminosity estimate is independent of the binary mass and takes into account that 3 — 10% of 
the total mass M of the binary is radiated over a time interval of ~ 100 M [38]. 
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situations - radiation inextricably linked to the past history that produced the inspiral. 
This is demonstrated in plots of extracted radiation from current simulations, where 
observers a distance R ext from the binary only see a flat radiation profile for the first 
t ~ -R cxt of evolution time. 

The post-Newtonian (PN) initial data developed in this and preceding papers [561 
[571 [T] is an attempt to address this shortcoming in initial data by incorporating to 
leading PN accuracy the gravitational-wave content of a physical binary inspiral. The 
previous paper in this series, (57], demonstrated the evolution behaviour of the data 
evolved with no numerical conditioning of the constraints. In summary, we noted that 
the extracted radiation agreed with physical expectations from the very start of the 
simulation (though the burst of junk radiation associated with puncture evolutions was 
not completely removed). We also encountered several related weaknesses in the new 
evolved data, including: 

• very high eccentricity (~ 10%) in puncture trajectories, until around 100M before 
merger; 

• (possibly related) an extremely slow stabilization of pre-merger horizon masses 
during evolution; 

• large constraint violations when all the gravitational-wave terms are included. 

In this paper, we present hybrid BY-PN data with partial constraint conditioning, aimed 
at resolving these issues. 

In Section [2j we summarize the standard conformally flat Bowen-York puncture 
initial data. In Section [3j we review the theoretical work that led to the PN data, and 
discuss the encoding of the past history of the binary in Section HI In Section [5J we 
present the new hybrid initial data, and demonstrate its improved numerical behaviour. 
The numerical evolutions and comparison for various initial data are discussed in 
Section |6j We conclude in Section [7J with a discussion of future development and 
application of this data. 



2. BY Initial Data 



Most groups use the puncture prescription of Brandt and Briigmann [55] with Bowen- 
York extrinsic curvature [51] to initialize the numerical fields for a puncture evolution. 
In this prescription, the three- metric 7y is conformally flat: 

7y = ip^dij , (1) 

and the conformal factor must satisfy the Hamiltonian constraint 

A^ + l^K^- 7 = 0, (2) 

o 

where the conformal Bowen-York extrinsic curvature already satisfies the momentum 
constraint for holes with arbitrary momentum and spin. Brandt and Briigmann's insight 
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to find the solution to equation was to factor out the divergent parts of ifi , leaving 
a well-behaved, simply-connected sheet on which to solve their modified constraint: 



where BY refers hereafter to Bowen-York puncture data, r a = |x — x a |, and m a and 
x a denote the "bare" mass and location of each hole, respectively. In this prescription, 
only a single elliptic equation for u has to be solved, and u will be regular everywhere on 
the grid. Despite the great success in simplifying the form of the constraint equations, 
this prescription does not make use of any information to accurately describe the past 
evolution of the binary black holes. Conformal flatness, for example, prevents the 
astrophysically expected gravitational radiation from being included in the initial data. 

3. PN Initial Data in ADM-TT Gauge 

Post-Newtonian (PN) techniques are considered to accurately represent an astrophysical 
system in the limit of slow motion/far- apart black holes. See Blanchet [58] and 
Schafer |5H] for a review of the PN approach. 

In the canonical Hamiltonian formulation with point-like sources in general 
relativity, Ohta et al [Stll EH E2] investigated a class of coordinate systems in 
which the metric tensor becomes Minkowskian at spatial infinity, and developed the 
2PN order calculation. In particular, they found a compact way of expressing a 
compatible coordinate condition with the aid of the ADM formalism. By applying 
a transverse-traceless decomposition of the three-metric and its conjugate momentum, 
they were able to encode the full information about the dynamics of the canonical 
fields and particle variables in a reduced Hamiltonian, which is obtained by solving the 
constraint equations only. To solve the constraints, we follow the steps of Schafer [63] . 
Historically, the 2PN order calculation was completed by Damour and Schafer [64] . 
Shortly afterward, Schafer [631 ES] included the 2.5PN radiation-reaction terms in this 
Hamiltonian approach. 

Tichy et al [56] adapted the 2.5PN ADM-TT results to puncture initial data for 
numerical relativity. We start then from the PN expression for the spatial metric, which 
differs from conformally flatness by a radiative term hJ T : 



^BY = 1 + E + U > 



(3) 



7™ = V4n<% + h% 



(4) 



where the PN conformal factor is given by 



# N = i + £^ + o(e 6 ), 

a=l Zr a 



(5) 



where e = 1/c is the PN order parameter and 




(6) 
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Here, p a represents the linear momentum of each hole, and r 12 is the separation between 
the holes. 

The accompanying extrinsic curvature K™ is related to the three- metric 's (trace- 
free) conjugate momentum 7Tp N : 



& = I7J- e 3 4 + e 5 4)+0( e *). (7) 



(See the complete equation (17) in [56].) Explicit expressions for tt%, and 7r^ can be 
found in [HE] ; The leading-PN-order term, 7rV 3 \ , coincides with the standard Bowen-York 
linear-momentum contribution for a nonspinning binary. 

Near each particle, the spatial metric can be approximated by: 

7™~(l + ^) 4 ^ + 0(l/r3), (8) 

which is just the Schwarzschild metric in isotropic coordinates. For r a — > 0, the 
coordinate singularity is approached. This represents the inner asymptotically flat end 
of the Schwarzschild metric in isotropic coordinates, or the puncture representation of 
the Schwarzschild solution. This shows that if the metric form is kept as in equation 01]), 
then there is a black hole centred on each particle. 

The traceless-transverse part of the spatial three-metric, on the other hand, can 
be constructed by imposing an outgoing wave condition. We can rewrite the evolution 
equation for hJ 3 T as: 

A? 



h TT = -D-l5 TTkl 

ij ret ij 



(9) 



where 5^ TiJ is the TT-projection operator as defined in [63] . Schafer [53] suggested a 
"near zone" approximation for h]j, by splitting the retarded inverse d'Alembertian in 
OH]) with an inverse Laplacian: 

hl?= -[A^ + M-A- 1 )]^"!...] 

= /l JT(NZ) + ^(remainder) +0(e5) _ (1Q) 

The explicit form of this near-zone (NZ) approximation was given by Jaranowski and 
Schafer [66J. A few years later, Kelly et al [57] completed the picture for nonspinning 
black holes by determining the "remainder" TT term, /^t (remainder) ^ 2PN order. 
The structure of the remainder term divides into three segments, according to time 
of evaluation: 

, TT (remainder) _ , TT (present) , TT (retarded) , TT (interval) -.n 
tl^j — Itij + llij + llij ■ { 11 ) 

For each field point where hj^ is to be evaluated, the "present" term is evaluated using 
the particle positions and momenta at t = 0. The "retarded" term is evaluated using 
positions and momenta at the retarded time of each source particle relative to the field 
point. The "interval" term is an integral over the particles' paths from the retarded 
time to the present. 



Hybrid black-hole binary initial data 



6 



The present-time piece almost completely cancels the near-zone solution. The 
kinetic terms (i.e., those involving particle momenta) cancel exactly, while potential 
terms (those involving particle-pair separations) are strongly suppressed. The retarded- 
time piece reduces to the quadrupole solution for a nonspinning binary as r 12 /r —> 
where r denotes the field distance. The interval piece is too difficult to do in generality: 
we must integrate numerically. 

The ADM-TT data introduced here has several attractive properties, which we 
briefly restate here. The three-metric and extrinsic curvature expressions are easily 
found. Unlike harmonic coordinates, ADM-TT has no logarithmic divergences. For 
a single BH, the data reduces to the Schwarzschild metric in isotropic coordinates. 
Up to 1.5PN order, the data coincides with the (unsolved) puncture approach - 
conformally flat, with Bowen-York extrinsic curvature. The trace of the extrinsic 
curvature vanishes up to 3PN order. The Hamiltonian constraint decouples from the 
momentum constraints. 

4. Encoding the Past History of the Binary 

To evaluate this initial data for a given separation at t = 0, we must know not only 
the particle momenta at t — 0, but also the position and momentum of each particle 
over the past history of the binary, at least from the retarded time relative to the most 
distant grid point in the numerical domain. That is, the larger the numerical domain, 
the further back in time we must look for the information needed to fill in the hj^ fields 
at these spatially distant points. 

To supply this position and momentum information, we evolve the Hamiltonian 
equations of motion, using a standard Taylor PN Hamiltonian and flux function 
evaluated at 3PN and 3.5PN order, respectively. The set of equations which we use 
in this paper has been given in [67J. This evolution is started at a suitably large 
separation (ri2 ~ 40M is easily sufficient for current purposes), and the separation 
and orbital phase are saved for later interpolation. We then apply a shift in time and 
orbital frequency to match the initial conditions of the fully numerical evolution. 

This shifted trajectory data is used directly to evaluate the present-time fields. We 
also use it to evaluate the two retarded times (one for each particle) at each field point 
by a nonlinear Newton-search algorithm. Then we can evaluate the retarded-time fields, 
and numerically integrate the interval terms. 

5. Hybrid Initial Data 

The similarities of the wavy ADM-TT leading-order terms and the Bowen-York puncture 
data may be the key to further constraining the initial data. We use the rationale that 
the largest Hamiltonian violation contributions are due to the PN conformal factor, ipp^, 
since hj^ is relatively much smaller for field points closer to the punctures. We propose 
then a hybrid initial data prescription in which we first solve the Hamiltonian constraint 
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for the traditional Bowen-York puncture initial data and then rescale and superpose the 
higher-order PN terms to this solution. More specifically, the Bowen-York extrinsic 
curvature is used as usual to source the Hamiltonian constraint for a conformally flat 
spatial metric: 

A^by + ^y^by = 0. (12) 

The solution follows then from the puncture trick and is given by 

1 .mi m 2 , „ \ 

V>BY = 1 + ^ — + — )+U, 13 

2 r\ r 2 

where u, the regular part of the solution, is determined by the numerical solution of the 
resulting Poisson-like elliptic equation. The functional similarity of this solution to the 
one provided by the PN expansion, equation (j5J), is quite striking. The singular part, 
or the Brill-Lindquist part [HE] if you wish, has the same functional form in each case, 
differing only in how the bare masses of particles and black holes are treated, and the 
missing regular contribution u in the PN conformal factor. In this new approach the 
ADM canonical quantities are then rescaled by the Bowen-York solution, i/)by, instead 
of the PN conformal factor ^pn : 

Hi = </>BY<% + C , («) 

*" = ^(«NS>+« , *8>). < 15 > 

where we note that corresponds exactly to the spinless Bowen-York contribution. 
The transverse-traceless part of the metric, fo£ T , and the higher-order PN contribution 
to the conjugate momentum, irlL, are then superposed with this puncture data. 



5.1. Constraint violations 



Looking into the violation of the constraint equations, Fig. [TJ it becomes clear that the 
PN metric becomes inaccurate close to the punctures, at least when compared to the 
traditional Bowen-York data. The long-dash-dotted brown curve there represents the 
constraint violation for the PN metric in the ADM-TT gauge, labeled in the figure and 
in all plots hereafter as "ADMTT PN". The purple vertical dotted lines indicate the 
apparent horizon location around one of the punctures; here the one located on the y- 
axis at y=— 4M. Focusing our attention on the plot on the left panel, the Hamiltonian 
constraint violation, we can observe a negative violation outside of the horizon for the 
PN data. This can be interpreted as, and has the effect of, an unphysical negative mass 
gravitating around each black hole. Unfortunately, the hybrid approach just introduced 
in the last section did not reduce this violation considerably, as the green dashed line in 
the plot shows. We decided then to investigate a way to improve the data approximation 
closer to the apparent horizons and that led to the second crucial element of the approach 
proposed here: the use of attenuation functions. 
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Figure 1. Constraint violations. These figures compare the violation of the constraint 
equations around a puncture located in the y-eotis at y=—4M. The (purple) dotted 
vertical lines indicate the location of the apparent horizon. The (red) solid line 
represents the violation for the traditional Bowen-York data. The long-dash-dotted 
(brown) lines describes the violation for the ADM-TT PN data, while the (green) 
dashed and the (blue) long-dashed lines correspond to the Hybrid and Hybrid with 
attenuation function, respectively. We keep the line styles consistently the same in all 
figures hereafter to facilitate comparison. The plot on the left panel represents the 
Hamiltonian constraint violation and the one on the right represents the momentum 
constraint violation. 



5.2. Attenuation function 

Since the PN approximation of the metric and extrinsic curvature is inaccurate for the 
inner or internal zone, i.e., the region around each black hole in binary systems where 
r i <§C r i2 an d r 2 r 12 , we need to consider a different approach to construct the initial 
data around that region. Johnson-McDaniel, Yunes, Tichy and Owen [69J introduced 
the black hole perturbation calculation by Detweiler [70] to be used as an approximation 
for the inner zone. 

Here, we treat the inner-zone metric by a simpler approach. We discard the PN 
contribution around each black hole by applying an attenuation function defined as 
follows (see Fig. [2] for its functional form along an interval on the y-axis): 



tanh In 



a t 2 mi 
ri 2 



\ at mi 

2\ - 1 / 



1 



tanh In 



r 2 



at m 2 



1 



a t 2 m 2 2 ' 
r 2 2 



(16) 



where at is a constant parametei|§|. Our hybrid initial data approach with the use of an 



§ Note that when we apply this attenuation function to , the correction at large distances is at 
higher PN order than we treat in this paper. At large distances, we a Taylor expansion with respect 
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Figure 2. The attenuation function along the y axis for a binary with mi = = M/2 
at yi = — £/2 = 4M, and at = 1. 



attenuation function then becomes: 




(18) 



(19) 



The effect of this procedure is a considerable reduction of the Hamiltonian constraint 
violation, as the blue long-dashed curve indicates on the left panel of Fig. [IJ 

While the Hamiltonian constraint violation is improved, and we are going to discuss 
its effect on the data evolution later on, the momentum constraint violation seems to 
be unaffected by these two key ingredients of our procedure. Both the Hybrid and 
the Hybrid data with attenuation seem to affect the momentum constraint residuals 
negligibly when compared to the PN data. We are also going to discuss later how 
we think this may be affecting the data evolution. But before discussing any further 
properties of this data, let us say a few words on how the numerical experiments were 
prepared and performed. This is the topic of the next section. 

6. Numerical Evolution 

The black-hole binary studied in these experiments had a mass ratio of q = 1, was 
nonspinning, with punctures initially located on the y-axis at ±4M. In this framework 
we use the puncture approach [55] along with the TwoPunctures thorn [71] to 
calculate the Bowen-York initial data. We evolved these black-hole-binary data-sets 
using the LazEv [72] implementation of the moving puncture formalism j3j H] with the 
conformal factor W = yfx = exp(— 20) suggested by [10] as a dynamical variable. For 
the runs presented here we use centred, eighth-order finite differencing in space [73] and 
an RK4 time integrator (note that we do not upwind the advection terms). 

to rrii/ri (i.e., a PN expansion) yields 



Therefore, the corrections due to the attenuation function are at 2PN higher order than the leading 
order of hJ, T . 




(17) 
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We obtain accurate, convergent waveforms and horizon parameters by evolving this 
system in conjunction with a modified 1 + log lapse and a modified Gamma-driver shift 
condition [74J |3], and an initial lapse a(t = 0) = 2/(1 + ?/> BL ). The lapse a and shift (3 a 
are evolved with 



where r\ is constant. 

We use the Cactus/Einstein Toolkit code [751 176] to provide the parallel 
infrastructure and the Carpet [77] mesh refinement driver to provide a 'moving boxes' 
style mesh refinement. In this approach refined grids of fixed size are arranged about 
the coordinate centers of both holes. The Carpet code then moves these fine grids about 
the computational domain by following the trajectories of the two black holes. 

We use the AHFinderDirect thorn [78] to locate apparent horizons. We measure 
the magnitude of the horizon spin 5* using the Isolated Horizon algorithm detailed 
in [79]. Once we have the horizon spin, we can calculate the horizon mass Mah via the 
Christodoulou formula [80] : 



where M irr = J A/(16n) is the irreducible mass and A is the surface area of the horizon. 



For the runs in this paper, we used a grid hierarchy spanning nine levels of 
refinement. We ran three different sets of grid hierarchy resolutions in order to determine 
the waveform accuracy, shown in Fig. [3j The low-, medium- and high-resolution 
hierarchies have mesh spacings of h c = {5M, AM, 3.2M} for the coarsest grid and of 
hf = {M/51.2, M/64, M/80} for the finest ones, respectively. The mesh spacings in the 
wave-extraction region are h cxt = {1.25M, 1M, 0.8M}. We apply reflection-boundary 
symmetry at the z = plane and 7r symmetry (that is, particle-exchange symmetry) at 
the x = plane. The outer boundary is located at 400M. 

The results comparing the different initial data sets, presented in the following 
sections, used the medium grid hierarchy only, i.e., the finest mesh spacing for this 
hierarchy is hf = M/64, resulting in approximately 30 grid points across the initial 
apparent horizon. 

6.1. Apparent Horizon Mass 

The first result from the evolution of the data introduced in the previous section that 
we would like to discuss concerns the apparent horizon mass. It was reported earlier [1] 
for the PN data that the mass was not being conserved with time - some sort of mass 
"leak" was plaguing the data evolution. We confirm this problem in Fig. HI as we can 
observe the long time it takes for the PN irreducible mass to relax towards a constant 
value. 



(d t -P % di)a = -2a K, 



(20) 




(21) 
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Figure 3. Phase accuracy. This figure shows the difference in phase for the (2, 2) mode 
of the radiative Weyl scalar ip4 at different resolutions. We compare the difference of 
the medium and low resolutions with that of the high and medium ones. We see that 
the phase error is below the threshold of 0.05 at MCI = 0.2 (vertical dotted green line) 
suggested by the NRAR collaboration [81) . 



The first step we took when trying to tackle this issue was to introduce the 
attenuation function for the PN data. As the purple dotted curve in the figure reveals, 
this approach alone was not successful. Only when we introduced the hybrid approach to 
the binary black-hole initial data, did we see an improvement in this unwanted feature. 
Now the mass conservation seems much more accurate, left only the presence of some 
bumps in the curve. 

We believe that the presence of these bumps in the data are related to unphysical 
negative mass, resulting from the initial Hamiltonian constraint violation. Since these 
bumps appear approximately every half orbital period, or roughly when the punctures 
cross the y-axis, it may be reasonable to assume that the initial Hamiltonian violations 
hangs in the domain around their initial locations or in the holes' "wakes" , occasionally 
being absorbed by the passing hole. The attenuation function seems to reduce this effect 
slightly, however only when associated with the hybrid data. 

6.2. Eccentricity 

Several methods have been proposed to estimate the eccentricity of numerically evolved 
binary inspirals; for recent examples, see [39| |82| [83] . We estimate the eccentricity of 
all four different data sets by first x 2 -fitting a polynomial function, r^it), to the binary 
coordinate separation, r NR (t), as a function of time. We subtract then the numerical 
data from the fitted function and normalize it by the same fitted function: 



The interpolating polynomial must be monotonically decreasing with no extrema in 
the fitting interval. We use the highest possible polynomial order satisfying these 




(22) 
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Figure 4. The behaviour of each hole's apparent horizon mass A/ah for Bowen- 
York, ADM-TT PN, ADM-TT PN with attenuation [(purple) dotted line], Hybrid and 
Hybrid with attenuation initial data. The mass "leak" was reduced considerably with 
the hybrid initial data. 

constraints. The assumption behind this procedure is a secular quasi-adiabatic shrinking 
of the orbital separation. The amplitude of the sinusoidal part extracted out of this 
orbital decay provides then an estimate of the orbital eccentricity. 

Specifically for these runs, where the initial separation is ro = 8M, we fit a second- 
order polynomial curve to the numerical data. We were careful to choose a time to stop 
fitting much earlier than the coalescence time. Usually we would eliminate the initial 
gauge effects from the data, and start fitting somewhere around t = 40M or later, as was 
the case for both instances of Hybrid data presented in Fig. [5j However this procedure 
turned out to be very sensitive to the initial fit separation for the Bowen-York and PN 
data. We decided then to include this initial transient in order to establish only an 
upper bound estimate for the eccentricity. We believe that this is not necessary for 
binaries initially farther apart. 

Amplitude readings from the eccentricity estimator plot as well as from a bare 
inspection of the coordinate separation as a function of time, Fig. |5l indicate that the 
Bowen-York data has an eccentricity of the order of 10~ 2 . Clearly the PN eccentricity 
is around 5 times this value. The Hybrid data alone, on the other hand, did not have 
a considerable impact on the eccentricity, possibly due to the fact that the momentum 
constraint violations were not improved as well. We could minimize a bit the high 
eccentricity effects by introducing the attenuation function, but this problem remains 
open. 
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Figure 5. Coordinate separation (left panel) and its eccentricity estimator in equation 
(right panel) versus time for Bowen-York, ADM-TT PN, Hybrid and Hybrid with 
attenuation initial data. 



6. 3. Waveforms 

We illustrate in Fig. O the presence of the realistic wave content in the hybrid initial 
data. The figure shows the z = slice of the data for the real part of the Weyl scalar 
^4 multiplied by the areal radius -Rareai- We show the traditional Bowen-York data for 
comparison in the right panel of the same figure. 

Like the wavy ADM-TT PN data, the hybrid data significantly reduces the 
amplitude of the spurious radiation when compared to the traditional Bowen-York data 
as is demonstrated in Fig. [7] for the real and imaginary parts of the (2, 2) mode of 
the Weyl scalar ijj^. In Fig. |8] we compare the amplitude (left panel) and phase (right 
panel) of this mode. We can observe a clear reduction in the amplitude and a dramatic 
improvement in the phase for all cases compared to Bowen-York data. Note that there 
is a decrease in amplitude for the Hybrid data, suggesting the binary is cast into a 
larger orbital separation. The use of attenuation function greatly reduces this effect, 
however. For the (4, 4) mode we only achieve a visible reduction of the junk radiation 
in the amplitude with the Hybrid with attenuation data, Fig. [9J which suggests that we 
need to improve the analytical approximation of the ADM canonical quantities in the 
inner zone if we want to accurately describe black-hole tidal effects. The disturbances 
in phase still seem greatly reduced. 

6.4. Final State 

Finally, we would like to comment on another advantage of the Hybrid data over the 
PN data. As Fig. [10] shows, the final state of the Hybrid data is much closer to the 
Bowen-York data than the PN one. Additionally, the post-merger horizon mass shows 
a significantly larger drift over time for the PN data than for the Hybrid data. 
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Figure 6. Comparison between the wavy pattern at t = present in i? area i7£e (^4) 
for the hybrid initial data (left) and the traditional Bowen-York data (right). 
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Figure 7. Real and imaginary parts of the (2, 2) mode of the Weyl scalar ip4 resulting 
from the evolution of Bowen-York, ADM-TT PN, Hybrid and Hybrid with attenuation 
initial data. The extraction radius is i? ox t = 90A-/. 



7. Conclusion 

Motivated by the benefits that the PN initial data could bring in bridging the fields of 
numerical and analytical relativity, we have introduced a hybrid approach to the initial 
data evaluation. We were able to significantly reduce the unphysical horizon mass loss, 
as well as reduce the eccentricity present in the PN data in ADM-TT coordinates, both 
undesired features previously reported in pQ. The method consists of solving for the 
traditional Bowen-York puncture data conformal factor and using this conformal factor 
to rescale the ADM canonical quantities as they appear in the PN data. The high 
eccentricity still present in the data was reduced by the use of an attenuation function 
in the black-hole inner zone, effectively decreasing the contributions of the higher-order 
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Figure 8. Comparison between the -04 (2, 2) mode amplitude (left) and phase (right) 
for the Bowen-York, ADM-TT PN, Hybrid and Hybrid with attenuation initial data. 
The extraction radius is -R oxt = 90M. 




Figure 9. Comparison between the -04 (4, 4) mode amplitude (left) and phase (right) 
for the Bowen-York, ADM-TT PN, Hybrid and Hybrid with attenuation initial data. 
The extraction radius is i? xt = 90M. 



PN terms and the transverse-traceless part of the metric to the canonical conjugate 
momentum and three-metric, respectively. 

We plan to continue our studies with two lines of research. First, as a short- 
term project, we would like to reduce the eccentricity of the hybrid data even further 
by exploring different attenuation functions and by the iterative eccentricity reduction 
procedure described in Pfeiffer et al [S3]. We plan to use larger orbital separations to 
facilitate this modeling and consequently the eccentricity reduction. 

The second plan of research, a long-term one, is two-fold. On the analytic 
side, we would like to explore alternative analytic approximations to the metric and 
extrinsic curvature in the binary black-hole inner zones, where the point-particle 
approximation breaks down. Our plan is to match the PN metric to a perturbed, 
boosted Schwarzschild black hole, a la Yunes et al [85], with a trumpet topology in quasi- 
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Figure 10. The apparent horizon mass and the dimensionless spin of the final black 
hole resulting from the evolutions of the Bowen-York, ADM-TT PN, Hybrid and 
Hybrid with attenuation initial data. 



isotropic coordinates [861 EH EE]; this should provide an analytic solution with smaller 
constraint violations. On the numerical side, we will implement a generic constraint 
solver for conformally curved data. By adopting a constraint-satisfying conformally 
curved metric we hope to better encode the past history of the binary and provide an 
initial data free of junk radiation. 
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